Giant Vortex Lattice Deformations in Rapidly Rotating Bose-Einstein Condensates 
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We have performed numerical simulations of giant vortex structures in rapidly rotating Bose- 
Einstein condensates within the Gross-Pitaevskii formalism. We reproduce the qualitative features, 
such as oscillation of the giant vortex core area, formation of toroidal density hole, and the pre- 
cession of giant vortices, observed in the recent experiment [Engels et.al, Phys. Rev. Lett. 90, 
170405 (2003)]. We provide a mechanism which quantitatively explains the observed core oscillation 
phenomenon. We demonstrate the clear distinction between the mechanism of atom removal and 
a repulsive pinning potential in creating giant vortices. In addition, we have been able to simulate 
the transverse Tkachenko vortex lattice vibrations. 
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Quantum liquids exhibiting superfluidity are able to 
support quantized vorticity, as has been often demon- 
strated in superconductors and helium superfluids. Re- 
cently, magnetically trapped gaseous Bose-Einstein con- 
densates (BECs) have been shown to accommodate vor- 
tices Q and arrays of singly quantized vortices In 
particular, large vortex lattices have been produced by 
rotating the condensate with an anisotropically deformed 
trapping potential Q and by cooling the initially rotating 
thermal gas directly into the vortex lattice state 0. 

Generically, multiquantum vortices (those with a wind- 
ing number larger than unity) are energetically disfavored 
and arc not expected to persist if created in rotating con- 
densates. However, a number of methods have been stud- 
ied in order to overcome this vortex dissociation instabil- 
ity, including an external repulsive pinning potential 0. 
Another, much studied possibility is to use steeper than 
harmonic trapping potentials to allow rotation of the con- 
densate at frequencies exceeding the centrifugal instabil- 
ity hmit of the harmonic trapping 0, QiHSEl- Man- 
ifestly multiply quantized vortices can also be created 
topologically by deploying the spin degrees of free- 
dom of the condensate accessible in optical traps Il2l — -a 
method which has also been realized experimentally [13j . 

Most recently, giant vortex structures have been pro- 
duced in rapidly rota ting harmonically trapped Bose- 
Einstein condensates |l4| . In this experiment, a near 
resonant laser beam was shone through the condensate 
in order to produce a density hole encircled by a high 
number of vorticity. The mechanism for the removal of 
the atoms from the beam volume was suggested to be the 
recoil from spontaneously scattered photons. Using this 
technique, oscillations of the giant vortex core area and 
preccssing motion of the giant vortices were observed. 

In this paper we simulate the experiment, and com- 
pare the effect of atom removal to that of a conservative 
repulsive pinning potential in the formation and stability 
of giant vortex structures. Our simulations are in quali- 
tative agreement with the observed experimental results 
reported in Ref. Il4l and we provide a simple quantitative 



explanation for the oscillational behavior of the giant vor- 
tex core area. In addition, we have been able to simulate 
the excitation of the transverse vibrational Tkachenko 
modes in the vortex lattice [l^l ITtI by applying the 
method used in the recent experiment [18|. 

We model the condensate with the usual Gross- 
Pitaevskii equation 
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dip{r, t) 
dt 



£(r,t)V'(r, t) 



(1) 



for the condensate wavefunction V'(r, t) in a frame rotat- 
ing with an angular velocity fJ. The operator 
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V2 + 14,t(r,t)-f t)p-f2-L 
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is the Gross-Pitaevskii Hamiltonian, with g determin- 
ing the strength of the mean-field interactions, and L 
is the angular momentum operator. The external po- 
tential yoxt(r,t) includes a harmonic trap, 14rap(r) = 
mu}\{x'^ + ?/^)/2, and an additional term MasGr(r,i) cor- 
responding to the laser beam used to perturb the equilib- 
rium state in the experiment. We model leaser as a local- 
ized beam whose time dependence is typically a square 
pulse, and the spatial shape is given by a Gaussian func- 
tion of the form 

11,cam(A,a;o,yo,r-o) = ^e-4[(--o)'+(^-^o)'l/-o. (3) 

Here the amplitude A is either real or pure imaginary, 
respectively corresponding to a conservative repulsive 
pinning potential or a localized loss term. The results 
presented here are obtained, unless otherwise stated, 
within 2D geometry using g = 1000 ftwj^flg, where Sq = 
^Jh/2muji_. We have initially chosen to normalize the 
wavefunctions such that the chemical potential fi = (£) 
determines the initial particle number. 

We start with an equilibrated vortex lattice rotating at 
a frequency U, = 0.9 uj_\_. An atom- removing loss term, 
Vbcam(— lOOi /iijJx, 0, 0, 7 oq), is then applied at the cen- 
ter of the trap for the interval, t = [0, 0.05T], where 
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FIG. 1: Oscillation of the giant vortex core area. The initial 
equilibrated vortex lattice rotates at angular frequency Q — 
0.9 u}±, and has fj, = 12huj_i_. Condensate density is locally 
suppressed from the vicinity of the trap center as described 
in the text. In effect, the giant vortex core area begins to 
oscillate followed by the excitation of the overall breathing 
mode of the condensate. The snapshots (a)-(h) are taken at 
times t = (0, 5, 19, 33, 60, 70, 81, and 100) x 10"^r. The sides 
of each picture are 40 5o wide. 

T = 2tt/ljj±. The subsequent evolution of the conden- 
sate density is displayed in Figs, nia)-(h) in the labo- 
ratory frame of reference. The area of the hole begins 
to oscillate with a dominant Fourier frequency 3.0 uj± 
during t = [0.05, 5] T, as illustrated by the lower oscilla- 
tion pattern in the inset of Fig. |21 [19j. The contraction 
and expansion of the hole reflects the hexagonal lattice 
structure, as is most clearly seen in Fig. ^d). The in- 
dividual vortex phase singularities within the giant vor- 
tices remain separated at all times in contrast to gen- 
uine multiply quantized vortex states. The loss of atoms 
from the trap center also excites the collective breath- 
ing mode(s) of the condensate mainly at the universal 
value 2.0 ujj_ ll^l, which is also obtained from our cal- 
culations for the Bogoliubov modes for the equilibrium 
vortex lattice eigenstates. In addition, both the core and 
breathing mode oscillations contain a second order In- 
dependent breathing mode component at 2.4 ujj_, which 
provides coupling between those modes. 

The core area oscillation of the giant vortex may be 
understood in terms of the centrifugal potential asso- 
ciated with the azimuthal condensate flow around the 
equilibrium radius of the vortex. The condensate flow 
field around the boundary of the giant vortex may be 
approximated to be w = — -, where £ is the number of 
circulation quanta enclosed in the giant vortex of radius 
r. Expanding the centrifugal potential in a Taylor 

series around an equilibrium radius Rq of the giant vor- 
tex, we find the quadratic term in the radial displacement 
to have the effective harmonic frequency 

UJcore = 2V3e(^^^ LU^, (4) 

suggesting the vortex density ^?/i?o to be the dominant 
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FIG. 2: Giant vortex core area oscillation frequency as a func- 
tion of the radius of the giant vortex. Data from our simula- 
tions ( X ) , joined by a solid line for clarity, are plotted together 
with the predictions (o) from Eq. ||1J. The dashed line is the 
breathing mode frequency 2 huj± for the 2D system. The data 
is obtained using the condensate shown in Fig.jSJa) as the ini- 
tial condition. The inset is described in the text, and the lower 
curve is scaled by a factor of 5. 




FIG. 3: Effect of applying a repulsive potential in the center 
of the vortex lattice, as described in the text. Frames (a)-(e) 
are for times t = (11, 27, 53, 71, and 100) x 10"^r. The initial 
configuration is as for Fig. Q. The sides of each picture are 
40 5,0 wide. 



factor determining the variation of giant vortex core os- 
cillation frequencies. 

To test the validity of this formula, we plot in Fig. |3 
our simulation results for the core area oscillation fre- 
quency as a function of the giant vortex core radius. The 
crosses (joined by a solid line to guide the eye) are the os- 
cillation frequencies obtained from the interval between 
the two first minima in the core size. The circles are the 
predictions of the Eq. 10} , in which we use for the ra- 
dius Rq the average of the first maximum and minimum 
core radii of the giant vortex. The number of circulation 
quanta £ is taken to be the number of vortices initially 
residing within the radius ro of the beam used to cre- 
ate the giant vortex . For smaller values of Rq (larger 
£/Rq) the core oscillates rapidly, but with increasing core 
size, the oscillation frequency slows, and approaches the 
value of the breathing mode, as the size of the giant vor- 
tex approaches that of the condensate itself. We have 



3 





FIG. 4: Angular momentum (upper curve) and condensate 
atom number (lower curve), as functions of time, normal- 
ized to their initial values. A weak atom removing loss term, 
Vbeam(— O.li fiwx, 0, 0, 5 ao), is constantly applied in the cen- 
ter of the condensate rotating at the frequency Q = 0.95 ui±. 
The solid squares are plotted for comparison with the experi- 
mental result, see Ref.ll4i The subplot shows a typical shape 
of the condensate after it has been strongly depleted. 



also performed 3D simulations for ft = 0.95 u!± with pa- 
rameters corresponding to the JILA experiment [l^ , and 
obtained an initial core oscillation frequency of 2.9 uj± in 
excellent agreement with the experiment |22j . 

Equation (0J also suggests that smaller rotation fre- 
quencies lead to higher core oscillation frequencies, since 
the equilibrium radius Rq of the giant vortex with con- 
stant vorticity isproportional to the healing length ^ cx 
(w^ — ri^)^^/^ [3. Howeyer, in the hght of our simu- 
lations, a more probable reason for the experimentally 
observed increase in the core oscillation frequency with 
decreasing rotation frequency is that the consequent in- 
crease in the vortex lattice spacing yielded giant vortices 
with smaller number of circulation quanta in them (as- 
suming constant radius for the apphed atom removing 
laser beam). For sustained oscillation, it is essential that 
there exists an energy barrier associated with the single 
vortices within the giant vortex 're-entering' the conden- 
sate. This holds the circulation quanta within the giant 
vortex together during its oscillatory motion. Eventu- 
ally, the condensate pushes back into the giant vortex 
core, and as we see in Fig.^ the first six vortices are lost 
back into the condensate from the corners of the hexagon 
at a time between the frames (d) and (e). If, instead 
of removing the atoms from the trap center, we apply 
a repulsive potential Vbcam(100 ?ia;j^, 0, 0, 7 oq), then, as 
shown in Fig. |31 no oscillation occurs in the giant vor- 
tex core area. Instead, the individual vortex cores are 
initially repelled from the trap center together with the 




FIG. 5: Formation of a toroidal hole in the rotating con- 
densate. Localized off-center Gaussian beam suppressing the 
condensate density is kept still in the laboratory frame. As 
the vortex lattice rotates, a ring is cut in the condensate. Af- 
ter one trap cycle T the loss term is switched off, after which 
the ring persists in the condensate for several condensate ro- 
tation periods. Frames from left to right are taken at times 
t = (0.05, 1, and 5) T. 



condensate and depending on the amplitude of the beam 
used, shock-wave-likc wave fronts and fine interference 
fringes may be produced. After removing the pinning 
potential, the mean-field pressure quickly closes the cen- 
tral hole leaving the whole condensate in the breathing 
mode oscillation. 

These simulations confirm that the mechanism for cre- 
ation of giant vortices in Ref. depends in an essential 
way on the removal of atoms within the volume of the ap- 
plied laser, and is different from the pinning effect from 
a repulsive laser beam. The reason is that while the pin- 
ning potential conserves angular momentum L and atom 
number A'', the absorbing beam increases L/N. This can 
be understood by noting that for the rotating lattice, the 
angular momentum cx flR'^ within a circle of radius i?, so 
that removal of atoms within R leaves an annulus with 
higher L/N. This annulus can now begin to relax towards 
an equilibrium vortex lattice, but in the absence of any 
thermal cloud, it must conserve the value of L, since the 
potential is cylindrical. It is straightforward to show that 
in comparison to the original lattice, a new lattice with 
less atoms and higher L/N must have larger angular ve- 
locity O. Thus as the annulus spreads into the central 
hole, driven by mean field pressure and confining trap, 
it encounters an angular momentum barrier because the 
condensate can not speed up rotationally as fast as the 
hole is filling. The core area thus begins to undergo oscil- 
lation. For the pinning potential on the other hand, even 
though atoms are repelled from the central region, the 
condensate retains the appropriate L and N for a vortex 
lattice rotating at the angular velocity fl. Hence, when 
the repelling potential is removed, the core fills imme- 
diately with no oscillation. The distinction between the 
two mechanisms is further emphasized in Fig. 21 where we 
plot the time evolution of the angular momentum of the 
condensate and the number of particles in the conden- 
sate while a weak atom removing loss term is constantly 
applied. The ratio 10/35 between the lost fractions of 
angular momentum and particles reported in the experi- 
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FIG. 6: Fundamental transverse Tkachenko vortex lattice 
oscillation mode excited by a localized loss of atoms, (a) 
The initial vortex lattice configuration with Q — 0.95 (i;x and 
= 11 hu!± is disturbed by applying a weak atom-removing 
loss term, Vbeam(— 0.02i ftoix, 0, 0, 4 So), for one trap period T. 
Consequently, the vortex lattice is set into slow vibrational 
motion around the rotation axis. Figures (a)-(c) are at times 
t = (0, 3, and 9) T. 



mcnt is in fair agreement with the data of our numerical 
simulations (at the time when 10% of atoms have been 
lost). Changing the amplitude A of the loss in the sce- 
nario of Fig. (0J, merely rescales the time. Furthermore, 
the shape of the curves shown are only weakly dependent 
on g. The subplot displays how the remaining parts of 
the condensate tend to break into a number of blobs af- 
ter the condensate has been intensively suppressed for a 
long time. Notice, especially, how the individual vortex 
phase singularities (-I-) are distributed over the whole gi- 
ant vortex core area. We find in these simulations that 
the centermost region (the low density area), also con- 
tains a number of anti- vortices (— ), which may be due to 
numerical noise. Moreover, it is worth pointing out that, 
unlike an application of a pinning potential, the method 
of atom removal cannot lead to stable giant vortex config- 
urations as all the atoms in the condensate would even- 
tually be consumed by the atom removing beam. 

In Fig. |3Ia)-(c) we display three frames of the rotat- 
ing vortex lattice when a localised off-centered atom- 
removing loss term Vbeam(100 fejj,, 6 aq, 0, 4 cio) is ap- 
plied. After one trap rotation period. Fig. [Sfb), the 
beam, stationary in the lab frame, has gouged a ring- 
shaped hole in the condensate. Depending on the radius 
of the applied loss term this giant vortex ring may per- 
sist for several condensate rotation periods. In a similar 
manner, if the beam is kept on only for a short period of 
time, the resulting giant vortex begins to precess around 
the trap center along with the rest of the vortex lattice. 
Again, if a pinning potential is used instead of removing 
atoms, neither dark rings nor precessing giant vortices 
are formed since the suppressed density heals soon after 
the influence of the pinning potential vanishes. 

We have also performed simulations using various val- 
ues for the damping term 7 in a dissipative version of 
the Gross-Pitaevskii equation, which models the effect of 
a co-rotating thermal cloud [23, |2J| . Using the damped 
model does not change the qualitative results here pre- 
sented but the time scales change and the giant density 



perturbations tend to heal more quickly. 

Finally, we have studied the recently experimentally 
observed transverse Tkachenko oscillations in the 
vortex lattice 0, 0, 0| . By suppressing the conden- 
sate wavefunction in the trap center with a localized loss 
term, we are able to induce transverse oscillations in the 
vortex lattice structure, as demonstrated in Fig. 1^ [25|. 

In conclusion, we have simulated the formation of gi- 
ant vortex structures in rapidly rotating Bose-Einstein 
condensates, and provided a simple mechanism which 
quantitatively explains the observed core oscillations. We 
have shown the clear distinction between the mechanism 
of atom removal and a conservative repulsive pinning po- 
tential in creating giant vortices. In the JILA experiment 
the formation of the giant vortex lagged the duration of 
the laser pulse, while in our simulations such a lag occurs 
only when the repulsive pinning potential is used. A pos- 
sible explanation of the lag in the experiment is that the 
repulsion caused by the atoms leaving the condensate 
creates outward pressure in the condensate delaying the 
giant vortex core formation. 
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